Quasi two-dimensional superfluid Fermi gases 
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We study a quasi two-dimensional superfluid Fermi gas where the confinement in the third di- 
rection is due to a strong harmonic trapping. We investigate the behavior of such a system when 
the chemical potential is varied and find strong modifications of the superfluid properties due to 
the discrete harmonic oscillator states. We show that such quasi two-dimensional behavior can be 
created and observed with current experimental capabilities. 



Finite size effects often lead to pronounced quantum 
mechanical changes in the system behavior. Among other 
things, these effects are crucial in the physics of quan- 
tum dots 1], in studying how spontaneous emission is 
modified in cavities 0, H 01 1 an d in understanding the 
properties of thin superconducting films . In addition, 
dimensional cross-over effects have been actively stud- 
ied experimentally using ultracold bosons in combined 
magnetic_and_pptical potentials as well as in optical 
lattices 



The possibility of fermionic superfluidity in ultra- 
cold atomic gases |l0( motivated a substantial exper- 
imental effort into controlling and manipulating cold 
fermionic atoms. Very recently, these efforts were re- 
warded with a spectacular success. A series of ex- 
periments demonstrated Bose-Einstein condensation of 




nance to vary the atomic interaction strength observed 
strong indications of fermion pairing. This was soon fol- 
lowed by the direct observation of the energy gap 
This set of experiments have established fermionic pair- 
ing in atomic gases as an experimental fact, and provided 
a strong case for superfluidity, proven recently by the ob- 
servation of vortices • 

The dramatic progress with ultracold fermions, com- 
bined with the fact that there is no fundamental problem 
in changing the dimensionality of the system by using for 
example optical lattices |2l[. raises important questions 
about the role of dimensionality and finite size effects 
in cold fermionic gases. The purpose of this Letter is 
to elucidate how quantum size effects are manifested in 
a dilute quasi-two-dimensional Fermi gas. Such system 
can be created by confining atoms tightly along one di- 
rection by a harmonic trapping potential. Note that we 
do not assume a purely Gaussian density profile in the 
z-direction, an assumption often made in studies of quasi 
two-dimensional quantum degenerate gases. Our exten- 
sion beyond the ground state Gaussian profile reveals dis- 
crete features in observable quantities, and is expected to 
be useful in studying how a true long range order at non- 
zero temperature is established as one crosses over from a 
purely two-dimensional into a three-dimensional system. 



We assume a cloud of fermionic atoms confined by a 
potential V(z) — muj 2 z 2 /2 and an atomic density n(z) 
which only depends on the z-coordinate. The cloud con- 
sists of equal amount of fermions on two different internal 
states denoted by | and |. Let us first, for orientation, 
consider the behavior of a non-interacting system at zero 
temperature: The first atom inserted into the system will 
fill the lowest energy state which corresponds to a station- 
ary atom in the xy-plane and whose axial wave function 
is the lowest harmonic oscillator state <fio(z) of the axial 
potential. This atom then has the energy Tiuj z /2. As we 
add more atoms, they fill the continuum of plane-wave 
states in the zy-plane with an axial wave function (j>o(z)- 
However, once the two-dimensional density of atoms is 
such that the next available plane-wave state has an en- 
ergy Tiuj z , the density of states changes abruptly and an 
extra atom has two choices. Either, it can occupy this 
high energy plane-wave state or it can occupy the zero 
plane- wave momentum state in the first excited harmonic 
oscillator state <fii(z). Similar doubling of choices occurs 
also with respect to higher harmonic oscillator states, 
when the continuum Fermi-energy matches the oscillator 
level energy. More quantitatively, the two-dimensional 
density (of both components combined) ri2D = J dzn(z) 
is related to the Fermi- momentum through k 2 F — 2irri2D- 
If the Fermi-energy Ep = Ti 2 kp/2ra is equated with 
the level spacing hto z we find a maximum density be- 
fore atoms start to occupy the first excited state of the 
harmonic oscillator, n c — l/(irl 2 ). (If not indicated oth- 
erwise, we use Ticu z as a unit of energy and l z = \Jh/mui z 
as a unit of length.) As the density exceeds this thresh- 
old, extra atoms have available a new state with the same 
total energy as the plane- wave states, but with a different 
axial state. This indicates that the growth rate of n2D 
with Ep is suddenly doubled from its earlier value of 
dri2D I 'dEp = m/(Trh 2 ). Therefore, dn^D I dEp for ideal 
fermions, at zero temperature, is a staircase with steps 
at Ep = Tilj z (n + 1/2) with the height m/iTrh 2 ). 

Let us now consider an attractive binary contact in- 
teraction with a coupling strength g between atoms in 
different internal states and allow for a non-zero tem- 
perature. Interaction between atoms is assumed to be 
sufficiently weak so that the Bardccn-Cooper-Schricffcr 
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(BCS) theory provides a reliable framework to study 
this many-body system also at finite temperature. If 
the Hamiltonian is expanded around the order param- 
eter A(z) = g{^)^{r)^ (r)), one finds the usual quadratic 
mean field Hamiltonian 
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This Hamiltonian is diagonalized using a Bogoliubov 
transformation $j(r) = J2c u c( r )^C,l + v l( v )^l i an d 

$j(r) = £^ — w c( r )^C,T + u c( r )^x l' wnere the quasipar- 
ticle amplitudes Uf(r) and Uf(r) are solutions to the 
Bogoliubov-de-Gennes (BdG) equations 



H A(r) 
A*(r) -H 



u c (r) y 4 V v c (r) 



where Hq = — h 2 V 2 /2m — fj, + Furthermore, the 

amplitudes are normalized J |w<;(r)| 2 + |w^(r)| 2 = 1. Self- 
consistency then imposes the well known gap equation 



A(z) = -.g]>> c (rK( r ) t 1 
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2n F (E c )} , (2) 



where n F (E) — 1/ (exp((3E) + 1) is the Fermi distri- 
bution. Finally the chemical potential is related to 
the atom density through the number equation n(r) = 
E c \u c (v)\ 2 n F (E c ) + K(r)J 2 (1 - n F (£ c )) . ^ 

In the quasi-two-dimensional system considered here, 
it is natural to expand the Bogoliubov quasiparticle am- 
plitudes in terms of the harmonic oscillator states 4> n {z) 
and radial plane waves ~ exp (ik ■ rj_). This amounts to 



u d r ) = EE:7^™( 2 ) e 



ik-r t 



(3) 



and the same expression for Wf(r), where A is the quan- 
tization area in the xy-plane. We include only the three 
lowest lying harmonic oscillator states and the solution 
to the BdG equations for the amplitudes u^ n , and k 
amounts to a diagonalization of the matrix 
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where £o,n,fc = huj z (n + 1/2) + 7?, k 2 /2m — /i, A„ 
j dzA(z)'|<Mz)| 2 , and A 02 = / dzA(z)<j)* 2 (z)Mz)- This 
matrix is almost block-diagonal with respect to differ- 
ent harmonic oscillator states. However, since 4'o( z ) an d 



cf>2 (z) are both symmetric, Ao 2 is non-zero and the simple 
block diagonality is broken. This coupling between n — 
and n — 2 channels is evident as an avoided crossing be- 
tween two of the three (positive) dispersion branches. 

For a weakly interacting system the restriction to just 
three harmonic oscillator states is expected to be suffi- 
cient, if the chemical potential is below the energy Thui z /2 
of the third excited state and if the temperature is small 
compared to hu> z . In our examples both these conditions 
are well satisfied. The coupling between n — and n = 2 
states influences, depending on /i, A(z = 0) by about 
10%. As a validity test, we included the n = 3 channel 
in the numerics and found a quantitative change of a few 
percent, while the qualitative behavior was unchanged. 

The eigenvectors w C,fc = 



' C C 



U l,k' V l,k> 



with eigenvalues E t 



*2,k> u 2,k 

can then be inserted into Eq. J5J) to find a self consistent 
solution for A(z). Since the Bogoliubov quasiparticle 
amplitudes are polynomials multiplied by the same 
exponential (in trap units) exp(— z 2 /2), it is conve- 
nient to expand the symmetric order parameter as 
A(z) = J2n=o a 2nZ 2n exp(— z 2 ). When the above expan- 
sion is used in Eq. (j2J) in conjunction with eigenvectors 
w^.fc, we find, by comparing terms with different powers 
of z, a set of three coupled self-consistency equations 
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(6) 



We solve this set of equations numerically. 

Replacing the sums over k with two-dimensional inte- 
grals gives rise to equations which are formally divergent. 
The ultraviolet divergence has its origin in approximating 
the interaction between atoms with a contact interaction. 
Most elegantly this divergence is removed by renormal- 
izing g to two-body scattering matrix, which amounts to 
subtracting the divergent part away from the integral. 
However, for computational reasons, we remove the di- 
vergence using a simple Gaussian energy cut-off. The 
cut-off energy is high enough so that our results are not 
sensitive to the cut-off procedure. 

The integrals are also infrared divergent and this di- 
vergence is due to the bound state appearing in the two- 
dimensional problem. This bound state has an energy 
eo /Tiujz ~ exp (— v27r£ z /lal V where a is the three dimen- 
sional scattering length 22]. This weakly bound state 
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causes a logarithmic energy dependence for the coupling 
strength g ~ In -1 (eo/e) which in turn results in non- 
separable gap equations. In our case this makes the com- 
putations exceedingly complicated. Furthermore, the 
bound-state energy also depends on the axial wave func- 
tion. This implies that the relevant bound state energy 
would also depend on the harmonic oscillator quantum 
number n and consequently one would have to deal with 
several different coupling strengths. However, in the BCS 
theory, the most interesting effects originate in the neigh- 
borhood of the Fermi surface. System behavior is there- 
fore largely insensitive to the details of the low energy 
behavior of the coupling strength. In fact, we found that 
the use of a sharp low energy cut-off is sufficient for the 
examples presented here. As long as the low energy cut- 
off was -C 1 , we could change it by an order of magnitude 
without affecting the results seriously. 

In Fig. ^ we show, for 6 Li atoms, an example of the or- 
der parameter in the center of the system A(z = 0) as a 
function of temperature and chemical potential. The gap 
A(0) increases with the chemical potential and shows a 
clear staircase structure especially around [i = 3htu z /2. 
The behavior around /i = 5Hlv z /2 is smoother. Also the 
critical temperature rises quite suddenly when the chem- 
ical potential is close to the harmonic oscillator levels. 
This increase reflects the abrupt increase of the density 
of states. While A(0) = ccq increases monotonically with 
chemical potential, this is not generally true for other co- 
efficients of A(z). Having solved the gap equations, the 
derivative dn2Djd[i can be easily computed. The stair- 
case structure one expects for an ideal Fermi gas is still 
present, but now the interactions have rounded the steps. 

In a homogeneous superfluid, the energy gap for the 
single particle excitations coincides with the order pa- 
rameter. In the inhomogeneous case, the order param- 
eter becomes position dependent and Andreev bound 
states (in-gap states) appear [2^. Note that in the RF- 
spectroscopy experiments so far, the final state was ini- 
tially empty, i.e. there is no Pauli blocking for transfer 
of particles from the lowest momentum states, unlike in 
superconductor-normal metal tunneling experiments. In 
a strongly interacting Fermi gas, however, the order pa- 
rameter can be as large as half of the Fermi energy and 
even the lowest momentum states are strongly affected 
by pairing. As a result, it turns out that the peak po- 
sition is of the same order of magnitude as the order 
parameter 0,0]- ^ n contrast, in the BCS limit, pairing 
takes place near the Fermi level and the transfer of par- 
ticles from the lowest momentum states leads to a very 
small shift of the peak position [25j], of the order A 2 /(2/i). 
Such small shifts are difficult to observe. However, in the 
near future it will be possible to trap stable mixtures of 
Fermions in three internal states of one atom (or, say, in 
one state of 40 K and two of 6 Li). Then, selected transi- 
tions can be Pauli blocked by preparing the Fermi level 
of the final state at will. This provides a new degree of 
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FIG. 1: Order parameter A(z) at z = as a function 
of dimensionless temperature T = HbT /TbU) z and chemical 
potential. We assumed a dimensionless coupling strength 
g — —0.25 here and in Fig. 2. 

control in the spectroscopy of the gas and the resolution 
at the BCS limit could be dramatically increased. In this 
way, also the Andreev state energies might be directly ob- 
served. An alternative way to resolve the lowest Andreev 
state energy could be via suppression of collective mode 
frequencies |26| . 

We calculated the RF-spectra for our system [25| and 
from that determined the location of the spectral peak 
as a function of chemical potential. This peak position 
and the lowest Andreev state energy together with the 
order parameter at z — are shown in Fig. [21 In the 
main figure, the peak position is calculated by assuming 
Pauli blocking of the final state, i.e. the final state chem- 
ical potential equals the initial state one. The spectral 
peak is closer to A(0) than to the lowest Andreev state. 
This trend becomes more clear as the chemical potential 
increases, meaning that the experimental signal in RF- 
spectroscopy from the lowest Andreev state becomes neg- 
ligible when several axial harmonic oscillator states are 
occupied. This might, however, be changed by the choice 
of the final state chemical potential. In the inset we show 
the peak position when the final state was initially empty. 
The spectral shift is now much smaller. However, both 
the main figure and the inset demonstrate that the steps 
in the order parameter are directly reflected in observable 
quantities. This behavior originates from the change in 
the density of states and demonstrates the many body 
nature of the pairing. 

Probing experimentally the regime we are interested 
in here requires a sufficiently tight axial confinement. In 
practice, the interesting two-dimensional density scale is 
ri2D ~ l/C 71 "^) an d when the axial profile of the cloud 
is that of the lowest harmonic oscillator state, this cor- 
responds to a density scale n^o ~ 1/(tvI^) at z = 0. If 
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FIG. 2: The lowest Andreev state single particle energy 
(dashed line), peak position of the RF-spectrum (dot-dashed 
line), and A(z = 0) (solid line) as a function of chemical po- 
tential at T = 0. The inset shows the peak position when the 
final state in the RF-spectroscopy is initially empty while in 
the main figure equal chemical potentials of the initial and 
final state were assumed. 

u> z = 2tt ■ 1000 Hz x s, where s is dimensionless, we get 
for 6 Li atoms n 3 o ~ s 3 ^ 2 x 1.5 ■ 10 n l/cm 3 . Assum- 
ing that the density in the real experiment is between 
10 13 1/cm 3 and 10 14 1/cm 3 , one would require axial trap- 
ping frequencies in the range of 10-100 kilohertz or more 
precisely s ~ 20 . . . 80. This is entirely realistic using, 
for example, a one dimensional optical lattice. Due to 
the larger atomic mass, for Potassium the required ax- 
ial trapping frequency is smaller by factor of about three. 
The features in Fig.[3]become of the order 1-10 kHz which 
can be resolved in an experiment. 

In an experiment the cloud is inhomogeneous also ra- 
dially. This inhomogeneouity is expected to broaden the 
staircases further. If the radial density profile varies on 
a length scale R, we expect a natural energy scale for 
inhomogeneity of ~ h 2 /2mR 2 and effects due to it be- 
come small if it is much smaller than the interaction en- 
ergy scale ~ |<?|7iw 2 = \d\fko z /l z . This implies a condition 
(R/l z ) 2 ^> (l z /\a\). Since the scattering length is tunable 
using Feshbach resonances this condition is not difficult 
to satisfy. Furthermore, many effects due to the radial 
density profile can be included, if necessary, using the 
local density approximation. 

In a purely two-dimensional system one expects 
a Berezinskii-Kosterlitz-Thouless transition associated 
with the proliferation of vortex-antivortex pairs. This im- 
plies non-trivial phase fluctuations of the order parameter 
A(z) and the absence of long range order (when T ^ 0) 
as required by the Coleman-Mermin-Wagner-Hohenberg 
theorem [i^. Here we have ignored such fluctuations 
with the implicit assumption that the modulus |A(z)| of 
the order parameter can be reliably calculated with the 
BCS theory. Understanding how the couplings between 



harmonic oscillator levels suppress phase fluctuations and 
establish true long range order would be an interesting 
future extension of the work. 

In summary, we have shown how finite size effects ap- 
pear and can be observed in a superfluid Fermi gas that is 
quasi two-dimensional via a tight harmonic confinement 
in one dimension. We introduced a theoretical approach 
which employs the first few oscillator states. Such an 
approach is sufficient to describe finite size effects, yet 
it is simple and transparent, which should be very use- 
ful in studies of phase coherence and fluctuations in di- 
mensional crossovers, of strong interactions, as well as of 
superfluid fermion dynamics in this system. 

We thank J. Kinnunen for helpful discussions. This 
work was supported by Academy of Finland (project 
number 106299). 
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